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We investigate tiie nontrivial interplay between geometry and temperature in the Casimir effect 
for the sphere-plate and cylinder-plate configurations. At low temperature, the thermal contribu- 
tion to the Casimir force is dominated by this interplay, implying that standard approximation 
techniques such as the PFA are inapplicable even in the limit of small surface separation. Thermal 
fluctuations on scales of the thermal wavelength lead to a delocalization of the thermal force density 
at low temperatures. As a consequence, the temperature dependence strongly differs from naive ex- 
pectations. Most prominently, thermal forces can develop non-monotonic behavior below a critical 
temperature. We perform a comprehensive study of such geothermal phenomena in these Casimir 
geometries, using analytical and numerical worldline techniques for Dirichlet scalar fluctuations. 



I. INTRODUCTION 

The Casimir effect inspiring many branches of 
physics 0, 01 J features a decisive geometry dependence: 
the fluctuation-induced interaction between test bodies 
or surfaces depends on their shape and orientation. This 
is because the Casimir effect arises from the fluctuation 
spectrum in presence of the surfaces relative to the vac- 
uum fluctuations. The spectral properties in turn are a 
direct consequence of the geometry. 

This geometry dependence becomes even more pro- 
nounced at finite temperature T: thermal fluctuations 
can predominantly be associated with a characteristic 
length scale, the thermal wavelength At ~ hc/{k-BT). 
Thermal fluctuations contribute to the Casimir force, 
whenever the scale set by the thermal wavelength is com- 
mensurate with a mode of the fluctuation spectrum as 
defined by the geometry. Therefore, thermal corrections 
to the zero-temperature Casimir effect generally cannot 
be described by universal additive terms or other simple 
recipes but require a careful analysis of the interplay be- 
tween geometry and temperature, as first anticipated in 

i- 

This "geothermal" interplay has first been verified in 
paradigmatic perpendicular-plates Q or general inclined- 
plates configurations Further evidence for the ex- 
perimentally relevant sphere-plate configuration has been 
provided recently in |7Hl(j|. Typical low-temperature de- 
pendencies in these open geometries obey power laws 
with characteristic exponents that are particular for the 
geometry. Most importantly, these power laws disagree 
with predictions from standard local approximation tech- 
niques such as the proximity force approximation (PFA) 
- even in the limit of vanishing surface separation. 
This is in contrast to zero-temperature forces which are 
often well described by the PFA in this limit jl^ . 

In this work, we perform a comprehensive study of the 
geometry-temperature interplay for the sphere-plate and 
cylinder-plate configuration. We study the Casimir forces 
induced by fluctuations of a scalar field obeying Dirichlet 
boundary conditions on the surfaces in order to explore 



the geothermal interplay in a most transparent fashion. 
Moreover, we use the worldline approach to the Casimir 
effect [l^ which on the one hand provides for a highly in- 
tuitive picture of the fluctuations, and on the other hand 
facilitates analytical as well as numerical computations 
from first principles [13413 ■ 

For instance, the failure of local or additive approxima- 
tion techniques can directly be inferred from the temper- 
ature dependence of the force density: the latter tends 
to delocalize for decreasing temperatures on scales of 
the thermal wavelength Q. Local approximation tech- 
niques may only be useful at finite temperature if the 
strict weak-coupling limit is taken [T8| . or in the high- 
temperature limit. 

In the present work, we analyze the thermal force den- 
sity distributions, compute thermal forces for a wide non- 
perturbative range of parameters, and determine asymp- 
totic limits. This facilitates a careful comparison with 
local approximation techniques, and, most importantly, 
yields new and unexpected results for the geometry de- 
pendence of thermal forces. For instance, the pure ther- 
mal force, i.e., the thermal contribution to the Casimir 
force, reveals a non-monotonic behavior below a criti- 
cal temperature for the sphere-plate and cylinder-plate 
case [13: the attractive thermal force can increase for 
increasing distances. This anomalous feature is triggered 
by a reweighting of relevant fiuctuations on the scale 
of the thermal wavelength - a phenomenon which be- 
comes transparent within the worldline picture of the 
Casimir effect. Whereas these non-monotonic features 
already occur for a simple Dirichlet scalar model, non- 
monotonicities can also arise from a competition between 
TE and TM modes of electromagnetic fluctuations in 
conflgurations with side walls (20l . l2l] | . 

While there are a number of impressive veriflcations 
of the zero-temperature Casimir force (2^ . a compari- 
son between theory and thermal force measurements suf- 
fers from the interplay between dielectric material prop- 
erties and finite temperature [l^, still being a subject 
of intense theoretical investigations (24| - |29t . In view of 
the geothermal interplay, we expect that the full resolu- 
tion of this issue requires the comprehensive treatment 
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of geometry, temperature and material properties, pos- 
sibly also including edge effects [30l - l33| . First results 
on the sphere-plate configuration using scattering theory 
and specific dielectric models demonstrate this nontrivial 
interplay 0,0,111. 

As a crucial ingredient for such an analysis, field- 
theoretical methods for Casimir phenomena have to be 
used that can deal with arbitrary Casimir geometries. 
In addition to the worldline methods [l^ [35l - l38| used in 
this work, a variety of approaches has been developed 
in recent years, such as a functional integral approach 
[sol - liH and scattering theory (i^ - fsoj . An extension of 
these methods to finite temperature is usually straight- 
forward and highly worthwhile in view of the geometry- 
temperature interplay. 

Our paper is organized as follows: after a brief ac- 
count of the worldline approach to the Casimir effect in 
Sect, nil the sphere-plate and cylinder-plate configura- 
tions are studied at zero temperature in Sect. IIIII In ad- 
dition to making contact with the literature, we perform 
the worldline computation directly for the force instead of 
the interaction energy. Section IIVI contains all our main 
results on the finite-temperature case. Our conclusions 
are summarized in Sect. [V] For reasons of comparison, 
the proximity-force approximation for the sphere-plate 
and cylinder-plate case is worked out in detail in ap- 
pendix 1^ In addition to explicit formulas which have 
not comprehensively appeared in the literature so far, we 
relate the PFA to an approximate treatment of the world- 
line path integral which helps to understand the differ- 
ences between the exact and approximate treatments. 

II. WORLDLINE APPROACH TO THE 
CASIMIR EFFECT 

We start with a short reminder of the worldline ap- 
proach to the Casimir effect for a massless Dirichlet scalar 
in 4 dimensions; for details, see [^, [l^ [s^l- For a con- 
figuration E consisting of two rigid objects with surfaces 
Si and E2, the worldline representation of the Casimir 
interaction energy in I? = 4 dimensional spacetime reads 

Ec = -^l^ % J d'^CM (es[x(r)]) . (1) 

The worldline functional 8s[x(t)] is 1 if the worldline 
x(r) intersects both objects S = Si U E2, and is zero 
otherwise. 

The expectation value in Eq. ([1]) is taken with respect 
to an ensemble of 3-dimensional closed worldlines with a 
common center of mass xcm and obeying a Gaufiian ve- 
locity distribution. For static Casimir configurations the 
time component cancels out at zero temperature. Equa- 
tion ID) has an intuitive interpretation: All worldlines in- 
tersecting both surfaces violate Dirichlet boundary con- 
ditions and are removed from the ensemble of allowed 
fiuctuations, contributing to the negative Casimir inter- 
action energy. During the T integration, the extent of a 



worldline is scaled by -s/T . Large propertimes 7" corre- 
spond to IR fluctuations, small T to UV fluctuations. 

Introducing finite temperature T = 1//3 by the Mat- 
subara formalism is equivalent to compactifying Eu- 
clidean time on the interval [0,/3]. The Casimir free en- 
ergy corresponding to Eq. ([T|) now becomes 
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For numerical purposes, it is convenient to remove the T 
dependence from the velocity distribution by the rescal- 
ing 

^(t):=^x(n) e-3/o^^''i^ = e-3J'o^'^*, (3) 
vT 

where 7 = d'y{t)/dt. Now, the 6 function reads more 
explicitly 

e[x] = e[xcM + \/r7W]- (4) 

The worldline integrals are evaluated numerically by 
Monte Carlo methods, i.e, the path integral is approx- 
imated by a sum over a finite ensemble of riL worldlines. 
Each worldline -^{t) is furthermore discretized by a finite 
set of N points per loop (ppl). To generate discretized 
worldlines with Gaufiian velocity distribution the v-loop 
algorithm was used in this work [isl. [sH. In the remain- 
der, we apply the worldline method to the sphere-plate 
and cylinder-plate Casimir configurations. 

III. CASIMIR EFFECT AT ZERO 
TEMPERATURE 

Let us first study the sphere-plate and cylinder-plate 
geometries at zero temperature. Here, we make contact 
with earlier results for the Casimir effect of a cylinder and 
sphere above a plate S [13, li^ - lii l. Moreover, we 
generalize the worldline method to directly compute the 
Casimir force instead of the energy which leads again to 
significant simplifications compared to previous energy 
calculations. The method will be generalized to finite 
temperature in the next section. 



A. Sphere above a plate 

We start the configuration of a sphere above a plate. 
The sphere of radius R is centered around the origin x = 
0. The infinitely extended plate lies in the 2 = — (a + i?) 
plane, where a is the minimal distance between both ob- 
jects, see Fig. [1] Since the configuration has a rotational 
symmetry with respect to the z axis, the three dimen- 
sional Xcm integration reduces to a two dimensional one. 



The Casimir energy ([T]) reads 

E^^-jQ^l :^yd^dzcMr(eE[x(r)]), (5) 

where we have switched to cyhndrical coordinates 
{r, zcm) with = x^m + Vcm- The 8E[x(r)] functional 
factorizes, 

es[x(r)] = es[xcM + \/r7]ep[xcM + Vfj]. (6) 

Here 0s and &p account for the intersection of a world- 
hne xcM + VT 7 with the sphere and the plate, respec- 
tively. Notice that 0s is independent of a, whereas 0p 
reads 

ep = 9{-{a + R + zcM + Vflz,^J), (7) 

where 7z„,i„ denotes the worldline's extremal extent into 
the negative z direction. 

As we are interested in calculating the Casimir force, 
Fc = —dEc/da, the derivative acting only on 0p pro- 
duces a d function which eliminates the z^m integral. 
The Casimir force thus simplifies to 

1 / f°° dT f°° \ 

F. = --(y^ -I d..es[s + vr7]). (8) 

Here, we have introduced 
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The transition from worldline calculations of the force 
does not only lead to technical simplifications. Also, the 
classification of relevant worldlines changes slightly: for 
the Casimir energy in Eq. ([T]) the worldlines are scaled 
by the propertinie VT with respect to their center of 
mass which is finally integrated over. For a given center 
of mass, all points on a worldline xcm + fiti) lie on 
rays originating from the center of mass. These rays are 
traced out by the T integral running from T = to 
T ~ oo. 

By contrast, the Casimir force in Eq. ([8|) results from 
worldlines which are attached to the point x on the plate. 
For a given point x, all points on a worldline x+ v'T7(ii) 
lie on rays which now originate from x. Again these rays 
are traced out by the T integral. 

Now, the plate is always touched by by construction 
for all values of T, the remaining problem being the de- 
tection of intersection events with the sphere. Adapting 
methods from [s^l, it is clear that only those points of a 
worldline lying on the rays intersecting the sphere even- 
tually pass through the sphere for some values of T. Let 
denote the set of points on such rays intersecting 
the sphere, with fc labeling these rays for a discretized 
worldline. Those values of propertime T for which this 
point lies exactly on the sphere can be obtained from the 
equation 

(X-f Vr7(ifc))' (10) 



3 




FIG. 1: Sketch of the sphere-plate configuration. The infinite 
plate (blue line) is in the z = —{R + a) plane. The sphere of 
radius R is in the origin. The center of mass xcm of the world- 
line is in (r, 0, —a — R — \fT 7z„i„). During the propertime 
integration the worldline always touches the plate, while all 
its points move on rays passing through (r, 0, —a — R). Only 
points lying inside the cone will pass through the sphere. 

Equation ([T0|) has two solutions, 

( si-^{t,) // x.7(tfc) y ixp-flA ' 
' [mtkW +V vi7(t.)iv mtkW ) ■ 

(11) 

For T € (7^~, Tj^) the point x+VT-yitk) lies inside the 
sphere. The point x can be viewed as a tip of a cone that 
wraps around the sphere with the opening angle 2a, with 
sin(Q!) = i?/|x|. The value of the square root in Eq. PT|) 
varies between zero and R. The square root is zero if 
the ray merely touches the sphere, and R if the ray lies 
on the cone's axis, i.e., if it coincides with the direction 
spanned by x. 

For a given r, the worldline intersects the sphere if 
the propertime T is in one of the intervals bounded by 
Eq. pT|) for all possible values of fc. Denoting these in- 
tervals by Tk ■= [T^ J '^^]- I'l^*^ total support of the 
propertime integral then is 

5(r)=|Jn. (12) 

k 

The r dependence of this support arises from the fact 
that the set of fc rays lying inside the cone depends on the 
position r where the worldline is attached to the plate. 
The Casimir force ^ now reads 

The most time-consuming part of the algorithm is the 
determination of S{r) if the distance between the sphere 
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and plate is small. To reduce the computational time, 
it is advisable to reduce the N points per worldline to 
the subset k < N points on the above mentioned rays 
intersecting the sphere. For a given r, all points on rays 
outside the cone can immediately be dropped. Further- 
more in the process of taking the r integral from zero 
to infinity, the opening angle of the cone shrinks. All 
points on rays which leave the cone through its upper 
half can then be dropped completely from the calcula- 
tion, as they will never enter the cone again. Only rays 
below the cone, i.e., between the cone and the plate, can 
enter the cone for larger values of r. With these opti- 
mizations and with one integral less, the computational 
time for Casimir force calculations is significantly reduced 
compared with those of the Casimir energies studied in 
previous worldline investigations. 

These simplification facilitate to extend the previously 
studied parameter range to even larger a/R ratios with 
higher statistics. 



B. Cylinder above a plate 

In many respects, the cylinder-plate configuration is "in 
between" the sphere-plate configuration and the classic 
parallel-plates case. This also holds for the experimental 
realization: the effort of keeping the cylinder parallel to 
the plate is less than it is the case for two parallel plates 
[5^ : for the sphere-plate case, this issue is simply absent. 
As a clear benefit, the force can, in principle, be made 
arbitrarily large, since it is proportional to the length of 
the cylinder. 

The geometry of the cylinder-plate configuration can 
be parameterized analogously to the preceding sphere- 
plate case: we consider the symmetry axis of a cylinder 
of an (infinite) length Ly and radius R to coincide with 
the y axis. The infinite plate lies in the z = —{R-\-a) 
plane, with a being the distance between the cylinder and 
the plate. 

The Casimir force can be obtained directly from 
Eq. ([T]), where use the fact that the 8s[x(t)] functional 
factorizes (cf. Eq. ([6])) 



0s[x(t)] — 6cyl[xcM [xcM + vTj]. (14) 

Here 6cyi a-i^d Op account for the intersection of a world- 
line xcM + VT 7 with the cylinder and the plate, respec- 
tively. Again, only Qp depends on a and is given in 
Eq. dll). 

The y integral in the Casimir energy ([1]) is now trivial 
due to translational symmetry. The Casimir force can 
then be obtained directly from Eq. ([8]) and reads 



where r 
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FIG. 2: Casimir Force of a sphere of radius R = 1 above 
an infinite plate vs. the distance a. The force is normalized 
to the zeroth-order PFA formula. We observe an excellent 
agreement with the exact asymptotic solutions for small a 
(44 ] and for large a [4^ up to a = 100. For larger a, the 
number of A'^ points per loop has to be increased far beyond 
N = 2 ■ 10^ used for this plot; otherwise, the sphere falls 
through the rough mesh provided by the insufficiently dis- 
cretized worldline, leading to a systematically underestimated 
force as is visible here for a > 100 (pink triangles). 



As in the case of the sphere, the worldlines x + -\/T7 
are attached to the plate at the point x. The only differ- 
ence is that the worldlines are now 2 dimensional - a fact 
which reduces the computational cost. Only those points 
of a worldline lying on the rays intersecting the cylinder 
pass through the latter for some values of T. The con- 
struction of the support of the T integral is identical to 
that for the sphere-plate case, such that the total Casimir 
force on the cylinder can be written as in Eq. ([Tc 
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dr I 



IS{r) 



r3 , 



(17) 



C. Zero-temperature results for the Casimir force 

It is instructive to compare our results not only with 
analytic estimates, but also with the much simpler prox- 
imity force approximation (PFA). The latter is used by 
default for the data analysis of geometry corrections in 
most experiments. It derives from a classical reasoning 
for generalizing the parallel-plate case; thus, deviations of 
the exact result from the PFA estimate also parameterize 
genuine geometry-induced quantum behavior. 

Roughly speaking, the PFA subdivides the surfaces 
into small surface elements, applies the parallel-plate 
force or energy law to pairs of surface elements and inte- 
grates the resulting force density. The PFA is inherently 
ambiguous as the measure for this final integration is not 
unique: possible alternatives are the surface measures of 
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FIG. 3: Casimir Force of a cylinder of radius R = I above 
an infinite plate vs. the distance a. The force is normalized 
to the zeroth-order PFA formula. We observe an excellent 
agreement with the exact asymptotic solutions for small a 
[44] and for large a Q up to o = 1000. 



one of the involved surfaces or any intermediate auxiliary 
surface. Later on, we will refer to the "sphere-based" or 
"plate-based" PFA as two generic options for the integra- 
tion measure. The PFA for the present configuration is 
discussed in detail in Appendix El 

The concept of the PFA can also be translated into the 
worldline picture: as an approximation to the ensemble 
of complicated multidimensional worldlines, we may re- 
duce the worldlines to one dimensional straight lines. The 
length of a line then corresponds to the average extent 
of a worldline into a certain relevant direction in a given 
geometry. This picture also explains the occurrence of 
deviations from the PFA as well as the sign of these de- 
viations in the Dirichlet case: due to the spatial extent of 
the worldlines, they generically intersect both boundaries 
for smaller values of T than simple straight lines. As 
small propertimes yield quantitatively larger contribu- 
tions, this property then results in a greater force. More 
precisely, the size squared of a worldline is proportional 
to the propertime parameter which is in the denomina- 
tor of the worldline formula, see Eq. ([1]). This explains 
why worldline results are typically underestimated by the 
PFA for small separations of the objects. For very small 
separations, the upper bound of the propertime integra- 
tion can effectively be set to infinity, whereas the lower 
bound is a measure for the first (proper-)time, when a 
worldline intersects both objects. 

The Casimir force for the sphere and cylinder is com- 
pared to the PFA estimates in Fig. [2] and [3] respectively. 
For similar comparisons for the Casimir energy, see [37| . 
We have normalized the force to the leading-order PFA, 
which is exact in the limit of vanishing separation a. The 
systematic study of the PFA, also from the worldline 
point of view, is summarized in Appendix [Xj We ob- 
serve that the normalized force obtained with worldline 



numerics does not lie inside the range spanned by the 
ambiguity of the PFA estimates. Most prominently, the 
sign of the deviations from the a — > limit is differ- 
ent in the Dirichlet scalar case, as can be understood in 
the worldline picture described above. These observa- 
tions have been frequently made in the literature before 
[l3,[3i,|4i,|4i,[5l. 

In order to obtain Fig. [5] we have used ensembles with 
up to riL = 1.6 • 10^ and 2 • 10^. At very small 

distances a the number of points per loop is not very 
important, since part of the systematic error is reduced 
by normalizing to the leading order result; thus, even 
= 5000 is sufficient for example for a = 0.0333 at a 
precision level of 0.1%. On the other hand at a = 100 
the number of points per loop used was 1.5 • 10^. For 
larger distances the number of points per loop has to be 
increased far beyond 2 • lO'', since, as shown in Fig. [2l 
Even such high resolution is not sufficient to resolve the 
small sphere for larger distances. 

Already anticipating our results for finite temperature, 
this observation gives us a rough estimate for the validity 
limits at small temperatures. Below, we observe that 
for a/i? <C 1 the maximum of the thermal contribution 
to the force density at low temperatures T < 1/R lies 
outside the sphere on scales r ~ 1/T. From the fact that 
ensembles with A^ = 1.5 • 10^ are reliable for those cases 
where the dominant contribution to the force density lies 
within r < 100, we conclude that temperatures above 
T > 0.01/i? are accessible also in the limit a — > 0. 

From an algorithmic point of view, the sphere-plate 
and cylinder-plate configurations differ with respect to 
computational efficiency also beyond the trivial dimen- 
sional factors: for a sphere at large separations, a large 
fraction of points of a worldline can be dropped right 
from the beginning, as they never see the sphere, i.e, 
they never lie on a ray inside the cone. The situation 
is different for a cylinder. Dealing with a two dimen- 
sional problem, we use two dimensional worldlines and 
the number of points per worldline which now have to 
lie in a wedge is higher than those lying in a cone for 
the sphere-plate case. Using comparable worldlines with 
a large number of points per loop, we thus expect the 
worldline numerics to break down at far larger distances 
a than in the case of a sphere. This is indeed the case as 
is visible in Fig. [31 

For the Fig. [3] we have also used ensembles with up to 
TT-L = 1-1 ■ 10^ and A^ = 2 • 10^. At a = 100, the num- 
ber of points per loop used was 3 • 10^, and increased to 
A^ = 1 • 10^ for a = 333 and up to A^ = 2 • 10^ for a = 1000. 
As expected, the required number of points per loop for 
a certain a is less here than in the case of a sphere. Even 
at such large separations as a = 1000, we observe an 
excellent agreement with (43j . The corresponding esti- 
mate for the validity limits at small temperatures then is 
T > 0.001/i?. 
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IV. CASIMIR EFFECT AT FINITE 
TEMPERATURE 

A. General considerations 

At finite temperature T — the free energy can 

be decomposed into its zero-temperature part Ec{0) and 
finite-temperature correction AEc{T), 



(18) 



The same relation liolds for the Casimir force Fc{T) = 

Fe(0) + AFe(r). 

Within the worldline representation of the free energy 
([2|) . the finite-temperature correction is purely driven by 
the worldlines with nonzero winding number n. Most im- 
portantly, the complicated geometry-dependent part of 
the calculation remains the same for zero or finite tem- 
perature. 

Let us first perform a general analysis of the thermal 
correction for a generic Casimir configuration following 
an argument given in We start from the assumption 
that the Casimir free energy can be expanded in terms 
of the dimensionless product aT, 



Ec{0) 



1 + ciaT + C2iaTf + csiaTf 



(19) 



No negative exponents should be present in Eq. (|T9| . 
since the thermal part of the energy disappears as T — > 0. 
Generically, the T ~ Casimir energy Ec{0) diverges for 
surfaces approaching contact a — !• 0. From Eq. ((T9)) . we 
would naively expect the same for the thermal correction. 
If, however, sufficiently many of the first c^'s in Eq. (fT9|) 
vanish, then the thermal part of the Casimir energy is 
well behaved and without any divergence for a — ^ 0. 

This turns out to be the case for two parallel plates 
(ci = C2 = 0, and i?c(0) ~ 1/a'^) and for inclined plates 
(ci = 0, and Ec{0) l/a^) Consequently, an ex- 

treme simplification arises: the low-temperature limit of 
the thermal correction can be obtained by first taking 
the formal limit a = 0. This was first observed in Q and 
then successfully applied in Q. 

In the following, we argue that there is no divergence 
in the local thermal force density in the limit a — ?> 
for general geometries. For a generic geometry, the a- 
divergent part can only arise from the regions of contact 
as a — > 0. The divergence for these regions at T = is 
due to the diverging propertime integral over 1/7"^+^/^ 
which is bounded from below by ^ a^. This is because 
for worldlines smaller than a the worldline functional is 
always zero. At finite temperature the divergence in the 
thermal correction for a is removed since one inte- 
grates now over exp(— n^/3^/4T)/T^^^/^, which is zero 
for every n > in the limit T — > 0. The only nonan- 
alyticity could arise from the infinite sum. That this is 
not the case can directly be verified: instead of integrat- 
ing over the support S, we integrate over T from zero to 
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FIG. 4: The (negative) thermal force density Eq. (f22|) for 
perpendicular plates (dashed blue line), cylinder above plate 
(dotted dashed line) and sphere above plate (solid red line) 
both of radius i? = 1 in the zero-distance Umit a — >■ for 
T — 1. The sphere-plate curve represents the radial density 
including the radial measure factor ~ 27rr. The thermal force 
densities of cylinder and perpendicular plates at r = are 
equal to the force density of two parallel plates, ~ 
0.1097. The thermal force density in the sphere-plate case has 
a maximum of « 27r x tt^ /90, where the factor 27r arises from 
the cylindrical measure. Note that a considerable fraction of 
the force density lies outside the sphere which only extends 
to r = 1. As the temperature drops, the maximum moves 
monotonously to the right. 



infinity, yielding 

roo exp y-^jS^ j 



E 

n=l 







7-I+-D/2 



dT = (2r)^r(D/2)C(D). (20) 



For finite temperature T > 0, Eq. ([20)) is a finite upper 
bound for the original local thermal force density. This 
procedure corresponds to substituting the critical regions 
of contact by broader (and infinitely extended) parallel 
plates, see 3- The thermal contribution is estimated 
from above by flattening the surfaces in the contact re- 
gion. The local thermal contribution to the Casimir force 
of the original configuration is clearly smaller than the fi- 
nite thermal contribution of parallel plates. As the latter 
does not lead to divergences for a — )• 0, there can also 
be no divergence for the general curved case arising from 
the contact regions. Of course, infinite geometries may 
still experience an infinite thermal force, as it is the case 
for two infinitely extended parallel plates, but the local 
thermal contribution to the force density will be finite. 
From a practical viewpoint, taking the limit a — > first 
simplifies the calculations considerably. 

Another important feature of low-temperature contri- 
butions to the Casimir effect is the spread of the ther- 
mal force density over regions of size 1/T even for 
very small separations a. This phenomenon has first 
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been demonstrated for the configuration of two perpen- 
dicular plates at a distance a (In this configu- 
ration, the sphere in Fig. [T] is replaced by a vertical 
semi-infinite plate extending along the positive z axis 
and an edge at z = 0.) The thermal force density 
A/c(r, T) = /c(r, T) — /c(r, 0) for this case as a function 
of the coordinate r on the infinite surface measuring the 
distance from the edge (i.e., the contact point at a = 0) 
can indeed be obtained analytically on the worldline from 
the thermal force, 

Here, Ai is a worldline parameter measuring the extent 
of half a unit worldline, i.e., the distance measured in x 
direction from the left end to the center of mass. It is 
clear from Fig. [T] that the lower bound in the T integral 
in Eq. (|2ip is given by r^/A^: this is the minimal scaling 
value for which the worldline intersects the semi-infinite 
vertical plate. From Eq. (PT|) . we read off the following 
force density: 

A/e(r,T) _ 7:^T^ 

Ly 90 




Analytic results for the thermal force can be obtained 
by rescaling the radial coordinate r — > Ai?' per worldline 
and using (Ai) ~ 7t/2. The thermal force between the 
perpendicular plates in the limit a — >■ upon integration 
then yields AFc{T) = -C(3) LyT^/iiT in agreement with 
i- 

The perpendicular-plates configuration is special as it 
features a scale invariance in the a ^ limit: Eq. ([22|) 
remains invariant under T — > Ta, r ^ r/a and A/c — >■ 
A/c/a'* for arbitrary a. As a consequence, knowing ([^^ 
for a single temperature value, say T = 1/a, is sufficient 
to infer its form for all other T. Equation H22p is shown 
for T ^ 1/R,R ^ 1 in Fig. [H For r < 1/T, the force 
density stays nearly constant, corresponding to the first 
term in (|22p. It rapidly approaches zero for r > 1/T. 
From this, we draw the important conclusion that the re- 
gion of constant force density in r direction can be made 
arbitrarily large by choosing sufficiently low T. 

Similar consequences arise for temperature effects in 
other geometries. We plot the thermal force densities for 
the sphere-plate and cylinder-plate configuration in Fig. 
m The thermal force density for a cylinder above a plate 
at a = has a shape similar to the one of two inclined 
plates, whereas the radial force density of a sphere above 
a plate exhibits a maximum due to the cylindric measure 
factor r, see Fig. [4l Although these force densities are not 
scale invariant due to the additional dimensionful scale 
R (sphere radius), its maximum nevertheless moves away 
from the sphere as the temperature drops. We conclude 



1 - 



o 




FIG. 5: Thermal Casimir force of a sphere above a plate 
in the limit a — )■ and for R = 1. The worldline result 
and the PFA predictions are normalized to the leading-order 
PFA, cf. Eq. ([26}. The leading-order PFA predicts a 
behavior of the thermal force for all T. On the other hand, 
for small T the plate-based PFA and the worldline result see 
a behavior, whereas the sphere-based PFA sees a ln(T) 
one. For large T ^ 1/R the behavior is for all curves. All 
predictions agree in the large T limit. For very small T, the 
worldline result runs into the blue area which spans the PFA 
predictions. 

that no local approximate tools such as the PFA will be 
able to predict the correct thermal force in particular 
at low temperatures. The fact that the force densities 
for sphere and cylinder are not scale invariant leads to 
different temperature behaviors for T < 1/R and T > 
1/R even in the limit a — > 0. 

B. Sphere above a plate 

Let us start with the expansion of the thermal force for 
a R and for small temperature T ^ 1/R. Following 
our general argument given above, no singularities in a 
appears in the limit a — >■ 0. Also, we expect that the 
thermal force decreases with decreasing R. This moti- 
vates an expansion of the thermal force with only posi- 
tive exponents for a and R. Assuming integer exponents, 
dimensional analysis permits 

AFc(T) ^ caRT^ + ciaT^ 

+ C2R^T* + csaRT'^ + O ((a/i?)^ {TRf) . 

(23) 

From our numerical results in the limit a — >■ 0, we ob- 
serve a behavior of the thermal force, see Fig. [5l We 
conclude that cq ~ is negligible with respect to C2 in 
the regime T > 0.01 where numerical data is available. 
In fact, we conjecture that vanishes identically cq = 0; if 
so, also ci vanishes, since the configuration would other- 
wise be more sensitive to temperatures at small a than at 



8 





4 








3.8 


:A : A-...r 
: : - ■■ A.. ., 


worldline i — i — i 
low T limit 




3.6 


r ' 




o" 

II 


3.4 






e 








fen 
< 

1 


3 
2.8 


i 





0.02 0.04 0.06 0.08 0.1 

T 



FIG. 6: Low-temperature behavior of -AFc(r)/r'* for the 
sphere-plate configuration in the limit a — >■ for R = 1. For 
0.03 < T < 0.1, we observe a linear behavior which can be 
fitted to AFc{a = 0, T) « -3.96 R^T'^ + 11.66 R^T'^. We have 
used 40 000 loops with 2 • 10^ ppl each. For T < 0.03, the 
number of points per loop used is not sufficient to resolve the 
sphere properly, inducing systematic errors (black triangles). 



a = 0. Our conjecture is supported by the following argu- 
ment based on scaling properties: the dimensionless ratio 
of the thermal correction and zero-temperature force has 
to be invariant under the rescaling 

a^a/a, R^R/a, T ^ aT. (24) 

The same holds for the ratio of AFc(a,T) at a = and 
the zero-temperature force at a ^ 0. For a <C i?, we can 
use the PFA for the zero-temperature force, which to 
leading order yields ^ R/a? . If cq 9^ 0, this leading ratio 
would be ^ co(ar)'^ which is invariant under the rescaling 
(p4| ; in addition, this ratio would be invariant under ((24|) 
with R fixed. If cq = 0, then this ratio is ^ C2Ra^T^ 
which is invariant only under the full transformation (|24ll . 

The result that for co 7^ the thermal correction would 
exhibit the same R dependence as the zero-temperature 
force for small distances a ^ i? is counterintuitive: 
whereas the radial force density in the small-distance 
limit at T is peaked right under the sphere near 
r ~ 0, the thermal correction arises from contributions 
at much larger r, cf. Fig. |31 As a simple estimate, we 
expect that the thermal correction is proportional to an 
effective area of the sphere, k, [a + R)2R + ttR'^/2, as 
seen by the worldlines. This estimate then is compatible 
with Co = ci = and C2/C3 « 1.8. 

The question arises why the PFA approximation yields 
a behavior despite the additional scale R. The reason 
is that R appears only in the combination r^/R in the 
force density, such that T aT, r — > r/a'^ leaves the 
force density invariant up to a multiplicative constant. 

Let us return to Eq. (|23p . For C2 and C3, we obtain 



numerically (see Figs. [S] and [Tj) 

C2 « -3.96(5), C3 « -2.7(2). (25) 

These numbers can be confirmed by the exact T-matrix 
representation [53|. Note that both coefficients have the 
same sign, implying that the absolute value of the ther- 
mal correction to the Casimir force increases with in- 
creasing a for sufficiently small a and T. This appar- 
ently anomalous behavior can be understood in geomet- 
ric terms within the worldline picture (l9| . 

The system has a critical temperature Tcr — 0.34(l)/i?: 
For T > Tcr, the thermal force decreases monotonically 
for increasing sphere-plate separation a in accordance 
with standard expectations. For smaller temperatures 
T < Tci-, the thermal force first increases for increasing 
separation, develops a maximum and then approaches 
zero as a — >■ 00. The peak position is shifted to larger a 
values for increasing thermal wavelength, i.e., decreasing 
temperature. In all cases, the force remains attractive, 
see Fig. [H As an example, room temperature T = 300K 
corresponds to the critical temperature for spheres of ra- 
dius R ~ 2.6/im. For larger spheres, room temperature 
is above the critical temperature such that the thermal 
force is monotonic. For smaller spheres, the thermal force 
is non-monotonic at room temperature. If, for instance, 
T = 70K and R = 1.6/im, the thermal force increases up 
to a ~ 9fim. 

The high-temperature limit T ^ 1/R agrees with the 
PFA prediction for a ^ and reads 

AFc(T^oo) =-^^r3. (26) 

In the limit a -> 0, the PFA yields Eq. ([26]) for all T. 
This is because geometrically the leading-order PFA cor- 
responds to approximating the sphere by a paraboloid, 
which is a scale- invariant configuration at a = 0. At 
finite a, the scale invariance is broken and a term ~ 
+R'K'^aT^ appears on the right-hand side of Eq. ([^ 
at low temperature in the leading-order PFA. By con- 
trast, we observe that the true a — > limit is charac- 
terized by a behavior for small T and behavior 
for large T. Also, the sign of the correction at finite a 
is different: the full worldline result predicts an increase 
whereas the PFA correction reduces the absolute value 
of the force, see Fig. [71 

It is interesting to compare our results to another PFA 
scheme beyond the leading-order PFA: the plate based 
PFA. This scheme is not scale invariant at a = 0, as the 
low-temperature limit for a <C i? is also quartic and given 

by 

AFP|A(a « i?, T ^ 0) = 7ri?2. (27) 

Equation ((27|) . in fact, corresponds to the thermal force 
density of two parallel plates integrated over the area of 
the region below the sphere, iiR^ . Numerically, the cor- 
responding worldline coefficient is more than ten times 
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FIG. 7: Plot of the a-dependent part of the (negative) ther- 
mal Casimir force for a sphere at a = 0.1 and i? = 1 in the 
low-temperature regime. In the range 0.03 < T < 0.08, the 
worldline data for -(AFc(0.1, T)- AFc(0, T)) is well approxi- 
mated by 0.268 - 1.57 (dashed blue curve), correspond- 
ing to a C3 coefficient in Eq. (|23p C3 ~ —2.7. Note that the 
absolute value of the Casimir force has increased with a. We 
have used 40 000 loops with 2 • 10® ppl each. 

larger than the PFA prefactor Tr^/gO « 0.345. In 
Eq. (|27p . the low-T behavior at finite a is exponentially 
suppressed, implying that the plate based PFA predic- 
tion for C3 is zero - which is again in contradiction with 
our worldline analysis. The formulae (pS)) and ^7} are 
derived in the Appendix. The thermal force at a = is 
shown together with the PFA predictions in Fig. [5l 

We now turn to the high-temperature limit, in a strict 
sense corresponding to T ^ 1/a and T ^ 1/R. The 
second requirement is automatically fulfilled in the small- 
distance limit a <C i?. Quantitatively, it turns out that 
the high-temperature regime is already approached for 
T > 1/a and T < 1/R. 

A special case arises for a — > 0, where the high- 
temperature limit agrees with the PFA prediction 
Eq. ([26)) in the leading order. For a > 0, the high- 
temperature limit is linear in T and the total force be- 
comes classical, i.e., independent of he. This behavior is 
rather universal being a simple consequence of dimen- 
sional reduction in high-temperature field theories, or 
equivalently, of the linear high-temperature asymptotics 
of bosonic thermal fluctuations [55l - l59| . In order to 
find the high-temperature limit, we perform the Poisson 
summation of the winding sum. The Poisson summation 
for an appropriate function / reads 

C50 00 

^ f{n/T)^V2^T ^ /(2^fcT), (28) 

n= — co fe= — cxD 

where / is the Fourier transform (including a l/-\/27r pref- 
actor) of /. Applying Eq. (pS)) to the winding sum, we 



FIG. 8: Thermal correction to the Casimir force of a sphere 
for various temperatures T and R — 1, normalized to the 
thermal force at a = 0. For sufficiently small temperatures, 
the absolute value of the thermal force correction AF{T) first 
increases with increasing a. For T < 0.05 the small a behavior 
is well described hy 1 + a{2.Q8RT'^ - R^15.7 T^)/ AF{a = 
0,r). This verifies the fit used in Fig. [7] with coefficient 
linear in ~ a in front of T^. From this prediction, we would 
expect the curves to be monotonically decreasing for T > 
2.7/15.77? « 0.17/7? in the leading order. Due to higher- 
order terms, the T = 0.2 curve still increases slightly first. 
The statistical errors of the worldline calculation are of the 
order of the thickness of the curves. 

obtain 

00 

= -1 + 2TVT^ + irVT^ exp {~T{2TTkTf) , 

k=l 

For finite a, the propertime integral is bounded from 
below and the last term is exponentially vanishing as 
T — >■ cxj. Evaluating the worldline integrals for the first 
two terms, we obtain 

Ai^e(a,T) = -F,{a) + TF,{a). (30) 
The evaluation of F{a) is analogous to Eq. (fT2| . 

where the support S{r) is the same as in the T ~ Q case, 
see Eqs. ^ and (HI]). 

The Casimir force remains attractive also for high tem- 
peratures. The function Fc(a), normalized to the PFA 
prediction 

(«) = (32) 
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FIG. 9: High-temperature coefficient Fc{a) for tlie sphere- 
plate configuration normalized to the corresponding PFA co- 
efficient. At large a, the normalized coefficient is conjectured 
to approach a constant. For small a, the behavior is well de- 
scribed by 1 -I- (0.14 ± 0.015) a. For a > 100, systematic errors 
similar to those of Fig. [5] set in (black triangles); at finite 
temperature, these errors are even more pronounced due to a 
softer propertime exponent 5/2. 
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FIG. 10: Thermal correction to the Casimir force for a sphere 
with R = 1 normalized to the zero temperature force for vari- 
ous temperatures and a < R. The worldline results (symbols) 
should be compared with the PFA estimate (lines). We ob- 
serve that this ratio of thermal to zero-temperature force is 
surprisingly well described by the PFA for a wide parameter 
range, especially in the high-temperature regime. This hap- 
pens because both Fc{a) and Fc{a) increase with respect to 
the PFA with roughly the same rate, see Fig. [5] and Fig. [21 



is shown in Fig. [51 The function a^Fc{a) is monotonically 
increasing on < a < 100 (similar to a^Fc{a)). At small 
a, we obtain 



FPF^(a) 



= 1-f (0.14±0.015)a. 



(33) 



In analogy to the zero-temperature force, we conjecture 
that also a^Fc{a) remains monotonically increasing and 
finally approaches a constant for a — > oo. A consequence 
of this conjecture is that the high-temperature limit then 
has a simple form, T 3> 1/a, without any relation to R. 
Indeed, demanding 



AFc(a,T) = -Fc(a) + TF,{a,T) < 



(34) 



for a fixed T, the limit a — > oo corresponds immediately 
to r ^ 1/a, since a^Fc{a) itself approaches a constant. 
Our numerical data shown in Fig. [51 is indeed compatible 
with this conjecture. However, the large-a limit is dif- 
ficult to assess due to the onset of systematic errors for 
a > 100. 

Comparing Fig. [2] and [9l we notice that the zero- 
temperature force Fc{a) and the high temperature coeffi- 
cient Fc{a) behave similarly. This is not surprising since 
Fc{a) in _D = 4 Minkowski space corresponds to Fc{a) 
in D = 3 Euclidean space due to dimensional reduction 
in the high-temperature limit. For finite a, the high- 
temperature limit is already well reached for T > l/2a. 
In the PFA approximation, the weaker thermal force at 
not too small temperatures is normalized by the weaker 
zero-temperature force, leading to an accidental cancel- 



lation, such that for T > l/2a 
Ai?(a)^_ 90C(3) 



Fc{a) 



■T 



aT « 3.49 aT, 



(35) 



independently of R. A comparison between the full 
worldline result and the PFA for the normalized force 
is shown in Fig. [TUl for various a and T. Since for small 
separations a < R, the PFA is a reasonable approxima- 
tion already at medium temperature l/2a > T > 1/2R, 
see Fig. [5l we observe that the ratio between thermal 
Casimir force and zero-temperature result is surprisingly 
well described by the PFA for quite a wide parameter 
range. We stress that the PFA is inapplicable for each 
quantity alone. 



C. Cylinder above a plate 

Analogous to the sphere-plate case, we start with the 
expansion of the thermal force at low temperature T and 
for a ^ i? as in Eq. (|24p . Again, we allow only for 
positive exponents for a and R. Even though ^/R terms 
appear in an a/i? <C 1 expansion at zero temperature, our 
numerical results at small finite temperatures, somewhat 
surprisingly, are consistent with an expansion of the type 



= C2RT^ 



(36) 



The potential leading-order terms cqVRT'^^^ and 
cw/aT''/'^ are expected to be zero similar to the sphere- 
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FIG. 11: Normalized thermal contribution to the Casimir 
force of a cylinder above a plate for i? = 1 in the limit a — > 0. 
The worldline result and the PFA predictions are normalized 
to the leading-order PFA Eq. ((3S}. The leading-order PFA 
predicts a t'^^^ behavior of the thermal force for all T. On 
the other hand, worldline numerical data is compatible with 
a T** behavior for small T ^ 1/R and a T"^ behavior for 
large T > 1/R. This is also observed in the plate-based PFA, 
whereas the cylinder based PFA goes as r''ln(T) for small 
T. All predictions agree in the high-temperature limit. The 
worldline result enters the blue area for small T into, which 
is the region spanned by the different PFA approximations. 



plate case, see above, since the configuration of a cylin- 
der above a plate is not invariant unter a ^ a/a and 
T aT. We have no evidence for a term ^ V aRT^, 
which would lead to a nonanalytic increase of the force. 
Thus at small temperatures, the powers of T are found 
to be integers in leading order. Similar to the sphere- 
plate case, we expect the low-temperature contributions 
to the thermal force to be proportional to the effective 
area ~ Ly{2R + a) seen by the distant worldlines. This 
results in the rough estimate C2/C3 « 2, which also im- 
plies that both coefficients have the same sign. 

In the limit a — 0, our data in the regime T > 0.01 is 
compatible with a behavior of the thermal force. For 
C2 and C3, we obtain (see Figs. [T2l and [T3|) 



C2 ~ -1.007(7), C3 w -0.41(4). 



(37) 



As in the case of the sphere, both coefficients have the 
same sign, i.e., the absolute value of the thermal Casimir 
force increases with increasing a for sufficiently small a 
and T < T^. For the critical temperature, we obtain 
Tcr w 0.31 (l)/i?. As in the case of a sphere, the ther- 
mal force decreases monotonically with increasing a for 



FIG. 12: Low-temperature behavior of -AFc(r)/r'' for a 
cylinder above a plate for 7? = 1 in the limit a — >■ 0. For 
0.02 < T < 0.06, we observe a linear behavior. In the 
range 0.03 < T < 0.05, our data can be fitted to the 
form AJ'c(a = Q,T)/Ly ^ ~1.0065RT^ + 3.163R'T^. For 
T < 0.02, systematic errors due to worldline discretization 
artifacts lead to a fast decrease of the data (black triangles). 
We have used 44 000 loops with 2 • 10® ppl each. 



The high-temperature limit T ^ 1/R agrees with the 
PFA prediction in the limit a — > as expected. 



AFc(r^oo) ^ 3C(l/2)C(7/2) /R 



T2 



-0.278v^r5. 

(38) 



As for the sphere, the PFA predicts the same force law 
([55)1 in the limit a ^ for all T. At finite a, the scale 
invariance is broken and a term ^ +0.185a\/i?T^/^ ap- 
pears on the right-hand side of Eq. (j38p in leading-order 
PFA at low temperature. By contrast, we observe dif- 
ferent power laws for different temperatures in the limit 
a -> 0: a behavior for small T and T*"/^ behavior for 
large T. Also, the sign of the finite-a correction of the 
full result is opposite to that of the PFA, see Fig. [131 a-U 
of which is reminiscent to the sphere-plate case. 

Incidentally, the beyond-leading-order PFA schemes 
reflect the correct behavior much better. We observe 
that the cylinder-based PFA turns out to be the better 
approximation (as for the sphere-based PFA in the pre- 
ceding section). This is opposite to the zero-temperature 
case. For the plate-based and cylinder-based PFA, we 
obtain 



90 



2i?w -0.219 i^T'' 



(39) 



T > Tcr; below the critical temperature, the thermal AfPFA(0,r^0) _ RT^ (3^^ + In (f^) + 180^(4)) 



force first increases up to a maximum and then decreases 
again approaching zero for a — >■ 00. The position of the 
maximum depends on T and increases with inverse tem- 
perature, see Fig. 1141 In both cases, however, the ther- 
mal force remains attractive. 



907r2 

(0.22 ln(i?r/27r) + 0.32)i?T''. 



(40) 



The plate-based result is equal to the thermal force of two 
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FIG. 13: Plot of the a-dependent part of the (negative) ther- 
mal contribution to the Casimir force for the cylinder-plate 
configuration at a = 0.1 and i? = 1 in the low-temperature 
regime. In the range 0.025 < T < 0.05, the worldline 
data for -(A_Fc(0.1, T) - AFc(0, T)) is well approximated by 
0.04125 r*- 0.335187 We thus conclude that the cs coeffi- 
cient in Eq. (|36|) is ~ —0.4125. Note that the absolute value of 
the Casimir force has increased with a. We have used 44 000 
loops with 2 • lO'' ppl each. 



parallel plates integrated over an area 2RLy. The plate- 
based coefficient is more than four times smaller than 
the worldline coefficient, whereas the leading coefficient 
of the cylinder-based formula becomes arbitrarily large 
as T 0. The formulae and (gOl) are derived in the 
appendix. The thermal contribution to the force in the 
limit a — >■ is shown together with the PFA predictions 
in Fig. HH 

Let us now investigate the high-temperature limit, 
which can be obtained by Poisson summation of the 
winding- number sum as in Eq. ([29|) . A special case arises 
in the limit a — > 0, where the high-temperature limit cor- 
responds to the PFA prediction Eq. ([55]) in leading order. 
For a > 0, the high-temperature limit is again linear in 
T and the total force "classical", i.e., independent of {he), 



AF,{a,T) ^ ^F,{a) +TF,{a) 
as in Eq. ([30| . we obtain 



87r3/2 



dr 



Sir] 



dr 

7-5/2 



(41) 



(42) 



where the support S{r) is the same as in the T = case, 
see Eq. P7|) . 

The Casimir force remains attractive also for high tem- 
peratures. The function Fc{a), normalized to the PFA 
prediction 



PFA 



3Vi?C(3) 
32V2a5/2 



(43) 



FIG. 14: Thermal correction to the Casimir force of a cylin- 
der above a plate for various temperatures T and R — 1 
normalized to the thermal correction at a = 0. For suffi- 
ciently small temperatures, the absolute value of the ther- 
mal force AF(T) first increases with increasing a. For T < 
0.05, the small-a behavior is well described by a(0.4125r* — 
3.35187 RT'^)/AF (a = 0,r). This verifies the fit found for 
Fig. [T3] with a coefficient linear in a in front of T^. From 
this form, we would expect the curves to be monotonically 
decreasing for T > 0.41/3.357? « 0.12/7? to leading order. 
Due to higher-order terms the T = 0.2 curve still increases 
slightly at the beginning. 



is shown in Fig. [TH The function a^/^Fda) is monoton- 
ically increasing for < a < 1000 and is reminiscent to 
a^^^Fc{a). At small a, we obtain 



At large a, we find using Eq 
temperature law (43| . 



1 + (0.125 ±0.017)a. (44) 
and the analytical zero- 



Jcja) 



1.46(2) 



,5/2 



(a + Rf ln(a + R) ' 



(45) 



We can compare our results with those of an analytical 
result [Hi in the limit R <^ H = R + a. The leading- 
order thermal contribution to the Casimir force in this 
computation based on scattering theory reads 



AF,{a,T)^LyT f 
Jo 



\n{qR) 



dq, (46) 



where the integrand has been approximated to leading 
order in \ar'^{qR). Here, st{q/T) is a 2'kT periodic saw- 
tooth function which in the range from to 2itT is given 
by -g/(27rT) 4- 1/2. The authors of [Hi have given a 
simple estimate of the integral for the limit i? <C 1/ 27rT 
by replacing ln((ji?) by ln(27ri?T) and carrying out the 
resulting integral. We compare our worldline results 
with Eq. (Uni) as well as with the simple estimate in Fig. 
[Tl]and[T71 
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FIG. 15: High-temperature coefficient Fc{a) normalized to 
the corresponding PFA coefficient for a cylinder above a plate. 
The large-a behavior is well described by Eq. (I45|) . At small 
a, we find a behavior ~ 1 + (0.125 ± 0.017)a. 



Here, we propose another estimate which is vahd for ar- 
bitrary T > l/(i?-|-a). In this case, the sawtooth function 
is approximately constant for q < l/(i?-|-a). We approx- 
imate the logarithm by inserting the value qq for which 
gexp(— 2g(_R + a)) is maximal: Qq = l/2(i? + a). In turn 
for T < 1/ (R+a), the logarithm can be approximated by 
insertion of the value go where (7cxp(— 2i3'(i?+a))st((7) has 
its first maximum: qo = ttT/2. We choose the first maxi- 
mum, as the integrand is oscillating for q > qo, such that 
cancellation can be expected to occur. However, choos- 
ing qo ^ T always leads to a regular T^/ ln(T) behavior 
for small T, whereas Eq. ([46|) changes sign at very small 
T, see Figs. [TH and [iTl We thus conclude that Eq. ^ 
is valid for not too small T. The thermal contribution to 
the Casimir force then reads 



AFe(a,r) 



-TLy d coth(27r(i? 



a)T) 



2Tr{R+a)T 



\n{qo) da 



8{R + a) 



(47) 



where qo = 2ttRT in the Emig et al. approximation 
[Hi, whereas qo = R/2[a + R) for T > 1/{R + a) and 
go = i?7rT/2 for T < I / {R + a) in the approximation 
proposed here. See Figs. [16] and [17] for the results at 
a = lOi? and a = lOOi? respectively. In the small T 
limit, Eq. (gT]) reads 



AFc(a,T) = 



2Tr^{a + R)T'^ 
' 451n(go) 



(48) 



Writing this as c{T,a, R)T'^, the coefficient c always 
disappears for go ^ as T — > 0. In our numerical world- 
line analysis, the systematic discretization errors lead to 
a vanishing of the corresponding coefficient as well, since 
the number of points per worldline becomes insufficient 
for resolution of the cylinder at very small T. For an 



FIG. 16: Thermal contribution to the Casimir force for a 
cylinder above a plate for a = 10 and R = 1 compared 
with the analytic result (|46p ("Emig et al." [i^ l and various 
approximations as discussed in the text. Here, we use the 
abbreviation H = R + a. Remarkably, our proposed estimates 
using go = ttRT/2 and go = 1/2//, cf. Eq. (|47)l . describe the 
actual behavior far better than the analytic result H46[l. which 
changes sign as T — >■ 0. Also the T > l/(7i-|-a) approximation 
using qo = 1/2H remains a reasonable estimate even for T < 
1/{R + a). For the worldline data, we have used 5000 loops 
with 2 • 10^ ppl and 7000 loops with 2 ■ 10*' ppl. 



increasing number of points per worldline, however, our 
data actually appears to point to a non-vanishing coef- 
ficient, see Figs. [TBI and [T71 In any case, we expect the 
leading-order multipole expansion which is behind the 
asymptotic result (|46p to break down at low tempera- 
tures due to the geothermal interplay. 

For large T on the other hand, Eq. (|47p becomes 



Ly {1 - T7r{a + R)) 
87rln(go)(i? + a)3 



(49) 



We observe that the negative of the T-independent part 
approaches the zero-temperature limit of the Casimir 
force for large a/R faster if we choose go — 1/{R + a) 
rather than go = l/2(_R+a). This choice of go ~ l/{R+a) 
then constitutes our second estimate for T > 1/{R + a). 

For not too small T, the analytic result and the various 
go approximations nicely agree with our worldline data, 
see Figs. [TBI and[T7l For higher temperature, the behavior 
becomes ~ T and the different results acquire different 
slopes which partly disagree for T — > oo. For a = 10, 
the analytic result becomes ~ 0.000424Lj,r, the go = 
l/{R + a) approximation yields ^ 0.000431ij^r, and the 
go = l/2(i? + a) approximation - 0.000334Lyr. The 
numerical worldline result is ^ 0.00041(3)Lj,T. 

At a = 100 and high T, the analytic result is 2.57 • 
10~^ LyT, the go = l/(i? -I- a) approximation yields 



jy-i ciiiu i-iic go = l/2(i?-|-a) approximation ■ 



2.66 •lO-^L.^T and the, 

''LyT. 

For large a + R, the temperature coefficient becomes 



2.31-10-%T. The worldline result is - 2.4(9) -10" 
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FIG. 17: Thermal contribution to the Casimir force for a 
cylinder above a plate for a = 100 and R = 1 compared 
with the analytic result (|46p ("Emig et al." j4^) and various 
approximations as discussed in the text. Even at such large 
separations, we observe that the analytic result (|46p becomes 
invalid and even changes sign as the temperature approaches 
zero. Incidentally our simple T > 1/{R + a) approximation 
using qo = 1/2H describes the actual behavior rather well 
also for smaller temperatures. For the worldline data, we 
have used 2500 loops with 6 ■ lO'^, 5000 loops with 3 • lO'' ppl 
and 14000 loops with 3 • 10^ ppl. 



0.125/(i? 



'ln(i? 



a) for both go = l/(i? + a) and 
l/2(i? + a). For the analytic result the corresponding 
prefactor is greater than 0.123 and may become 0.125 
for i7 — > oo. The corresponding worldline prefactor is 
0.116(2), see Eq. (|45l). 

Let us return to the high-temperature discussion, and 
finally remark that also in the case of a cyliirder the ther- 
mal Casimir force normalized to the zero temperature 
result is well described by the PFA for T > l/2a. Anal- 
ogously to Eq. (j35p . we conclude from the dimensional- 
reduction argument, that the ratio of thermal to zero- 
temperature force in the high-temperature limit T > 
l/2a is approximately 



AJ-c(a) ^,^ 72C(3) 
Fcia) 



■T 



aT K 2.79 aT. 



(50) 



Also at medium temperatures this ratio is surprisingly 
well described by the PFA, even better than in the case 
of a sphere, see Fig. [TOl The normalized thermal force is 
shown in Fig. [18] for various a and T. 



V. CONCLUSIONS 

In this work, we have analyzed the geometry- 
temperature interplay in the Casimir effect for the case 
of a sphere or a cylinder above a plate. Since finite- 
temperature contributions to the Casimir effect are in- 
duced by a thermal population of the fluctuation modes. 



FIG. 18: Thermal correction to the Casimir force for a cylin- 
der with 7? = 1 normalized to the zero-temperature force for 
various temperatures and a < R. The worldline results (sym- 
bols) are in a better agreement with the PFA estimates (lines) 
than in the case of a sphere, see Fig. 1101 The normalized PFA 
results agree with the worldline results for not too low T be- 
cause again both Fc{a) and Fc{a) increase with respect to the 
PFA with roughly the same rate, see Fig. [3] and Fig. 1151 



the geometry has a decisive influence on the thermal cor- 
rections as the mode spectrum follows directly from the 
geometry. A strong geometry-temperature interplay can 
generically be expected whenever the length scale set by 
the thermal wavelength is comparable to typical geome- 
try scales. 

Within our compreheirsive study of the Casimir effect 
induced by Dirichlet scalar fluctuations for the sphere- 
plate aird cyliirder-plate geometry, we observe several sig- 
natures of this geometry-temperature interplay: the ther- 
mal force density is delocalized at low temperatures. This 
is natural as only low-lying long-wavelength modes in the 
spectrum can be thermally excited at low T. As a conse- 
quence, the force density is spread over length scales set 
not only by the geometry scales but also by the thermal 
wavelength. This implies that local approximation tech- 
niques such as the PFA are generically inapplicable at 
low temperatures. Quantitatively, the low-temperature 
force follows a power law whereas the leading-order 
PFA correction predicts a behavior - a result which 
has often been used in the analysis of experimental data. 
Only for ratios of thermal to zero-temperature forces, we 
observe a potentially accidental agreemeirt with the PFA 
prediction for larger temperatures. Here, the errors in- 
troduced by the PFA for the aspect of geometry appear 
to cancel, whereas the thermal aspects might be included 
sufficiently accurately. 

Another signature of this geometry-temperature inter- 
play is the occurrence of a non-monotonic behavior of the 
thermal contribution to the Casimir force. Below a criti- 
cal temperature, this thermal force first grows for increas- 
ing distance and then approaches zero only for larger dis- 
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tances. This phenomenon is not related to a competition 
of polarization modes as in (20l . [2]] | , but exists already for 
the Dirichlet scalar case. The phenomenon can be under- 
stood within the worldline picture of the Casimir effect 
[l9| being triggered by a reweigliting of relevant fluctua- 
tions on the scale of the thermal wavelength. From this 
picture, it is clear that the phenomenon is not restricted 
to spheres or cylinders above a plate; we expect it to oc- 
cur for general compact or semi-compact objects in front 
of surfaces, as long as the lateral surface extension is sig- 
nificantly larger than the thermal wavelength. In fact, 
another consequence of the delocalized force density is 
that edge effects due to finite plates or surfaces will be 
larger for the thermal part than for the zero-temperature 
force. 

Our results have been derived for the case of a fluc- 
tuating scalar field obeying Dirichlet boundary condi- 
tions on the surfaces. For different fields or boundary 
conditions, the temperature dependence can significantly 
differ from the quantitative results found in this work. 
This is only natural as different boundary conditions can 
strongly modify the fluctuation spectrum. For instance, 
the thermal part of the free energy in the sphere-plate 
case exhibits different power laws for Dirichlet or Neu- 
mann boundary conditions in the low-temperature and 
small-distance limit For future realistic studies of 
thermal corrections, all aspects of geometry, tempera- 
ture, material properties, boundary conditions and edge 
effects will have to be taken into account simultaneously, 
as their mutual interplay inhibits a naive factorization of 
these phenomena. 



Appendix A: The proximity-force approximation 
(PFA) 
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FIG. 19: The effective heights predicted by the PFA for a = 0, 
R = 1 and T = compared with those of worldline numerics 
for a sphere and a cylinder, respectively, obtained from the 
T = force density. The PFA predictions lie in the blue area 
which is bounded by Eq. ()A4|) from below and by Eq. HA8|) 
from above. The effective heights seen by worldlines are well 
approximated by the leading order PFA for r not too large. 
We conclude that for small a/R the force is described best 
by the leading order order PFA, since for small a/R the force 
density is concentrated around r — 0. 



At finite temperature, the corresponding expression is 
AePP(/i) _ 



Cpp 



(A3) 



1 9or ^ 



CO 



th(2n7r/ir) -|- 2m:hTcsc\?{2mihT) 



The proximity force approximation is a scheme for esti- 
mating Casimir energies between two objects. In this ap- 
proach, the surfaces of the bodies are treated as a super- 
position of infinitesimal parallel plates, and the Casimir 
energy is approximated by 



Ql) dcr. 



(Al) 



Here, one integrates over an auxiliary surface S, which 
should be chosen appropriately. The quantity e^^(/i) de- 
notes the energy per unit area of two parallel plates at a 
distance h apart, which at zero temperature reads 



_PP 



(A2) 



where cpp — 7r^/1440 for the Dirichlet scalar case. 

As the PFA does not make any reference to boundary 
conditions, all the formulas in this appendix are analo- 
gously valid for the electromagnetic case; all force formu- 
las then have to be multiplied by a factor of two for the 
two polarization modes. 



The distance is conventionally measured along the nor- 
mal to E. The two extreme cases in which S coincides 
with one of the two bodies provides us with a region span- 
ning the inherently ambiguous estimates of the PFA. 

For a sphere at a distance a above a plate, we thus inte- 
grate either over the plate ('plate based' PFA), or over the 
sphere ('sphere based' PFA). The Casimir force is then 
obtained by taking the derivative of (jAip with respect 
to a. However, for the 'sphere based' PFA, dh/da ^ 1 
(see below). This implies that deriving the force estimate 
from the PFA of the energy in general is not the same as 
setting up the PFA directly for the force (the latter would 
correspond to a surface integral over the parallel-plates 
force per unit area). In this work, we use the derivation 
via the energy (jAl[) . 

The dependence of the PFA prediction on the choice of 
E disappears in the limit a — >■ at zero temperature to 
leading order. This result shall be called 'leading-order' 
PFA. It can also be obtained by expanding the surface 
of the sphere/ cylinder to second order from the point of 
minimal distance to the plate and then using the 'plate 
based' PFA for this expansion. 
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factor approaches 1 for small r, all functions h coincide 
in this limit. 

The PFA can also be developed within worldline for- 
malism. Calculating the Casimir force density for two 
parallel plates, we have to determine that value of prop- 
ertimes T for which one dimensional worldlines, attached 
to one of the plates, touch the other plate for the first 
time. This event is encoded in the lower bound of the 
proper time integral, whereas the upper bound is set to 
infinity. Thus, we obtain 



4,T 



327r2 



'(VA)%,= -cx 



r3 



-dr . (A9) 



FIG. 20: The effective heights predicted by the PFA for 
a = 10, R — 1 and T = compared with those of worldline 
numerics for a sphere and a cylinder, respectively, obtained 
from the T = force density. Note that the effective height for 
the cylinder is in a local maximum at r = 0. At greater sep- 
arations the heights seen by worldlines are on average lower 
than the PFA predictions resulting in greater Casimir force. 
Also at larger separations the leading order PFA reflects best 
the actual situation, while the 'sphere/cylinder based' PFA 
turns out to be the worst. 



The representation (jA9P is suitable for zero and low tem- 
peratures, whereas for high temperatures one should use 
in (|A9[) the Poisson resummed winding sum 1^^. We 
encounter cumulants of worldline extents A in low and 
high temperature limits which can be determined via the 
analytic expression (STj 

(X^) =D{D-l)T{D/2KiD). (AlO) 



The corresponding expressions for h read 

hpB = a + R- ^JW~^, (A4) 

n-SB = /iCB = TTT - n., (A5) 

cos(t^) 

^LO=« + ^. (A6) 

For /ipB, we integrate over [— i?, i?], for /ilo over all r 
and for /isB, ^CB over [— 7r/2,7r/2] with an appropriate 
measure. Note that right underneath the sphere/cylinder 
all h are equal to a. Demanding d^^ = dr for 6* ^ we can 
transform the integration over Q into an integration over 
r in a simple way by the substituting sin(0) — > r j R. The 
integral then goes from — i? to i?, and the corresponding 
h reads 

/isB = /icB = -R + R^^^. (A7) 

Also a measure factor resulting from Rd9 ~ R/\/R^ — 
and d/i/da have to be taken into account. At zero tem- 
perature, we can absorb these factors into the new effec- 
tive height 

With or without the prefactors ft-sB/CB is always greater 
than /ipB and /ilo ^-nd diverges for r — > i?. Since the 



Let us now point out the difference between the PFA 
and the worldline approach. In the PFA, we always use 
one dimensional worldlines to determine the distance, 
whereas the worldline dimension in the full formalism 
corresponds to the dimension of the geometry. To obtain 
the Casimir force for configurations containing one infi- 
nite plate in the worldline formalism, we integrate over 
this infinite plate as in the plate-based approach. How- 
ever, the integration does not stop at the end of the sec- 
ond body, which in the present case is a sphere or cylin- 
der. At arbitrary large distances, there are still worldlines 
which see the sphere/cylinder, i.e., we have to integrate 
to infinity. We therefore expect the leading-order PFA 
to refiect best the exact force laws. However, the proper- 
time support is not the same, and thus worldlines see an 
effective height different from the one of the leading-order 
PFA, see Fig. [19]and[20l 

The shape of the effective worldline height is roughly 
the same for zero and high temperatures. But at low tem- 
perature, the worldlines are reweighted. Only worldlines 
for large propertimes contribute considerably and thus 
worldlines at larger distances from the sphere become 
increasingly more important. Also their inner structure 
comes into play. Using the 'plate based' PFA, we ignore 
these effects and take into account only the region below 
the sphere/cylinder with the same function hpB ; hence, 
the result is expected to be too small. 

In the following, we apply Eq. (jA9P (multiplied by 
dh/da if necessary) to find the PFA expressions for the 
sphere and cylinder above an infinite plate, respectively. 
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1. Sphere above a plate 

a. Leading- order PFA 

For the sphere, the evaluation of the leading-order PFA 
results in an especially simple expression, 



d pPFA d /""^ 



r dr 



da 



da 



'o (a + r^/2RY 



(All) 



27ri? cpp 



da 



dh 



'*LO 



LO r, n PP/ \ 

27rite (a). 



Obviously, the relation F^q^ = 2ttR e^^ (a) remains valid 
also at finite temperature. We thus obtain 



rSHa,T = 0) = - 



271 R Cpp 



-K^R 

720a3 



(A12) 



At finite T and small aT {aT ^1/2), Eq. ([All yields 



AFLT(a,r) = - 



i?C(3) 



aR-K^ 
45 



(A13) 



For large aT (aT ^ 1/2), the expression ((X3| leads di- 
rectly to 



(a,r) 



TT^R 

720a3 

pPFA 



i?C(3) 



F--(a,0)+TFLT(«) 



(AM) 
(A15) 



Note that at a = the leading-order PFA predicts a 
behavior of the thermal force for all T. At finite a, 
the validity of the low-temperature limit is independent 
of R. With increasing a the absolute value of the PFA 
thermal force is always reduced, irrespective of T, quite 
the contrary to the full worldline results as discussed in 
the main text. 



Plate-based PFA 



Using Eq. (|A9p . we obtain 




Let us first analyze Eq. (|A16P for a — 0, 



AF^i^ 




(a-0,T) = - 
e^iT^R--! R^ 
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(A17) 
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Equation (jAl7[l distinguishes low- (T<Cl/i?) and high- 
temperature (T^l/i?) regimes. The low-temperature 
regime is already well approached for T ^ 1/2R. For 
higher T, the thermal force is in the high-temperature 
regime, T > l/2i?. 

At low temperatures, we have a behavior which is 
given by the first term in Eq. (|A17p . For higher T, this 
term is canceled by the term with the exponen- 
tial function, such that the leading behavior is given by 
the term. Then, expanding Eq. (|A17p . we get a 
contribution: 



AF^i^ia^O,T)^ 
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C(3)i? 
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Subtracting Eq. (Pl8| from Eq. (|XT7| and performing 
the Poisson resummation, we obtain the full T^\/2R 
behavior at a = 



ATpT 



(o,r) 



C{i)RT^ 



72 

C(3)r 

8R 



1440i?2 ■ 



(A19) 



Thus, the leading large-T behavior at a = is T''. 

Let us now consider the case a ^ 0. For a R and 
low temperature T ^ l/2(i? -I- a) « l/2i?, we have a 
behavior given by the first term in Eq. (|A16[) . The depen- 
dence on a is exponentially suppressed. This corresponds 
to the case of two parallel plates with an area of ttR^, 
where the dependence on a is suppressed exponentially 
as well. 

At medium temperature, 2{R + a) w 2R Z 1/T k.2a 
only the second and third term in Eq. (jA16|) are expo- 
nentially suppressed and can be neglected. The leading 
order can be found by expanding the remainder and con- 
sidering only the converging sums. 

To find the subleading terms, we again perform the 
Poisson resummation. For medium temperature 2R ^ 
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1/T ^ 2a and a <Si R, we then obtain 

a{a + 2R)t:'^T^ (a + R)C,{i)T^ 



(a,T) 



90 



(A20) 



C(3) 



72 8(a + i?) 1440(a + i?)2' 

The high-temperature hmit for 1/T ^ 2a can be per- 
formed irrespective of the actual value a/ R by summing 
up the whole Eq. (jA16|l . The result reads 



7r3i?2(3a + 2i?) i?^ C(3) 
1440a3(a + Rf 

(a,0) 



-T (A21) 



pPFA^ 
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^{a + RY 
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Note that Eq. \A2l\ reduces to Eq. (|A14|) for a 
it should. 



as 



For a larger than a ^ R, the following temperature 
behavior occurs. At low temperature 1/T ^ 2(i? + a), a 
T'^ behavior arises from the first term in Eq. ()A16p . At 
higher temperatures, the behavior becomes rapidly linear 
as given by Eq. (|A20p . being valid for 1/T ^ 2a. 



The plate-based force can be obtained in closed form 
from Eq. (|A3p also without using the worldline language: 
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Sphere-based PFA 



For the sphere-based PFA, the thermal Casimir force is given by 



R 



2 i-tt/2 



AF^i^ia, T) = -(—J^ s\n{9)h's^ia)d9 



''sb/-^^ n=l 



(A24) 



where /isB(a) is given by (jA5|) . For a = 0, we obtain the PFA approximation using the worldline language 
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where 7 is Euler's constant. The expansion in T does 
not terminate after a few terms, so we concentrate on 
the two leading coefficients. The coefficient in front of T'^ 
contains the worldline average (In A) . For an analytical 
expression, we note that InA = mlnA^/™. For large m, 
we get A^/'" 1, such that we can expand the logarithm, 



(InA) = ( lim m(Ai/"- 1)) 



7/2 + ln(27r), 



(A26) 



where we have used Eq. (|A10|) . Thus, the small-T limit 
of Eq. ((X25)) reads 



AFii^{0,T) 



j?2(2C^(4)-FC(4)(3 + 21n(j f:))) 

TT 

- R^C{5)T^- 



(A27) 

At a = 0, the PFA estimate | AFgPf^(a = 0, T)| fies above 
|A^]P|A(a = 0,T)|, see Fig. El For not too small T, the 
worldline result lies above both these PFA predictions, 
but due to the logarithm in the T'^ coefficient the sphere- 
based PFA becomes larger at smaller T, such that the 
worldline force enters the area spanned by the PFA pre- 
diction, see Fig. [51 
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The high-temperature hmit can be obtained by ex- 
panding Eq. ()A25|) about T = oo. The converging terms 
give the leading-order behavior. For the subleading or- 
ders, one has to perform the Poisson summation, how- 
ever, the integral involved is rather complicated and may 
still be inflicted with artificial convergence problems. The 
leading-order behavior for a = and large T reads 



AFii^O^T) = -^T^ + + OiT), (A28) 



and corresponds to the leading behavior of the plate- 
based hmit (|A19|) . 

Let us turn to the case of finite a. Expanding Eq. 
(jA24|) . we obtain the a dependent part of the thermal 
force, 
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The series has a form of (i?/a) ln(l -f- a/R), which we 
verified explicitly to 10th order. Assuming that this form 
holds to all orders, we get 



(a,r) 



(0,T)._1.. 



(A30) 



Note that the first two terms agree with the coeffi- 
cient of the plate-based formula (jA20p ; we also see that 
the absolute value of the thermal force decreases with in- 
creasing a. As Eq. (jA30P was obtained by interchanging 
summation and integration, we cannot expect Eq. (jASOp 
to describe the full a dependence for all a and T. Indeed 
at a fixed, the thermal correction AFgQ^{a,T) becomes 
~ T as T ^ oo, which is clearly not the case for Eq. 
(lASOl) . 

We can estimate the range of applicability of Eq. (jA30p 
as follows. At high temperature and a sa 0, all PEA 
estimates agree. For large T, the leading behavior is 



^ T^, see e.g. Eq. (|A28p . With increasing a the force 
is still attractive. Demanding AFgQ^{a,T) < we see 
that that Eq. (jA30|) leads to a positive thermal force for 
a > 1/T. On the other hand in the low-temperature 
regime, the a = contribution is given by Eq. (jA27p . 
Taking only the leading contribution into account and 
demanding AF^^^{a,T) < 0, we again obtain that the 
force becomes positive at a > 1/T. These rather rough 
estimates demonstrate that the validity range for a be- 
comes narrower with increasing temperature. For very 
small a, however the thermal correction is linear in a ir- 
respectively of T, whereas the dependence on a in the 
plate-based PEA is exponentially suppressed for small 
T< l/(i? + a). 

At large temperatures T > 1/a, we have the familiar 
situation 

AFii\a, T) = -Fii\a, 0) + TFii\a), (A31) 
where 



^3 (6a2 - iaR + 2R^) In (^) 



lAAOa^R 



240i?2 



and 



F, 



PFA 
SB 



(a) 



(i?-2a)C(3) ln(^jC(3) 



8a2 



4i? 



(A32) 



(A33) 



2. Cylinder above a plate 



a. Leading order PFA 



Unfortunately, a simple relation similar to F^q^ = 
27ri?e^^(a) does not hold any longer for the cylinder, 
such that the resulting formulas are not related to the 
known results of parallel plates and are rather compli- 
cated. For arbitrary a and T, we obtain 



J 



pPFA, 



a,T) 




4' 4' 



2 ' 4a2r2 J 
I 



F \ ^ ^-^ 2- 



4a2r2 




(A34) 



where 2F2 is the hypergeometric function in the standard 
notation. Eq. (jA34p does not distinguish between a < R 
and a > R, since the relevant parameter for different 
temperature regions is aT. For small aT {aT ^ 1/2), we 



can expand Eq. (jA34p . resulting in 



F^^^ia^T) _3Vi?C(7/2)C(l/2) 



y7/2 



15av^C(9/2)C(-l/2) 
4V2n 



(A35) 
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For large aT, the Poisson resummation of Eq. (|A34|) leads 
to 



32\/2a3 



768%/2a4 



(A36) 



which is, of course, -Fl^^{a, 0) + TFl^^{a). Note that 
at a = the leading-order PFA predicts a T^/^ behavior 
of the thermal force for all T. At finite a, the validity 
of the low-temperature limit is independent of R. With 
increasing a, the absolute value of the thermal force is 
always reduced, irrespective of T, quite the contrary to 
the full worldline results. 



and 



pPFA 
^PB 



(a) 



3(l + a)C(3) (3 + 4a + 2a2)C(3) 



16a5/2(2 + a)5/2 
3(1 + a)ArcTan 



Va(2+a) 



2(l + a)(2 + a)27r 
C(3) 



3a5/2(2 + a)5/27r 



(A39) 



b. Plate-based PFA 



Here, we give only the analytic expressions for special 
limits, since no general expression could be found in a 
closed form. At a = and T ^ ^/R, the thermal force 
can be found from the result of two parallel plates with 
an area oi A = 2RLy, 



In Eqs. ((A38|l and ((A39| . we set i? = 1; general expres- 
sions can be reconstructed by simple dimensional analy- 



sis. 



c. Cylinder-based PFA 



AF^E^{a = 0,T) = -2RLy-T\ 



(A37) 



As temperature rises, the T behavior changes from T"' 
to poj. T ■> 1/R and a = 0, the plate-based PFA 

agrees with the leading-order PFA, and the thermal force 
is given by the first term in Eq. (jA35|) . At low temper- 
atures and a <C i?, the dependence on a is exponentially 
suppressed, just as in the case of the plate-based PFA for 
the sphere. 

Finally, at finite a and T ^ 1/a, the force becomes 
classical -F^i^ia, 0) + TF^i^ia), with 



(15 -f 2a(2 + a)(ll + 3a(2 + a)))Tr'^ 



1440a3(l-^a)2(2-f a)3 
(5 + 4a(2 + a))7r3 



960aV2(2 + a)V2 

(5 + 4a(2 + a))7r2ArcTan 



(A38) 



480aV2(2-t-a)7/2 



For the cylinder-based PFA, the thermal Casimir force 
is given by 



AF 



PFA 
CB 



7r/2 



E 



dr. 



r3 



(A40) 



where hcB{a) is given by Eq. (|A5p . At a = 0, the thermal 
force can be found in closed form, 
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- 1260V2r 
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(A41) 



where /„ is the modified Bessel function of the first 
kind, pFq the generalized hypergeometric function, c = 
2RT^/nSf_imd x = nX/2RT. At small T, the expansion 
of Eq. (|A41|) leads to 

AfPg^(0,T) _ i^T^ (37r4 + 2^^ In (||^) + 180C^(4)) 



907r2 

6^ 



(A42) 



As temperature rises, the T behavior changes to T^/^^ 
For T 1/R and a = 0, the cylinder-based PFA agrees 
with the leading-order PFA and the thermal force is given 
by the first term in Eq. (|A35|) . For sufficiently small a 
the difference to the a = result reads 



AFEE\a,T)-AFSiHO,T) 



an 
'45 



T^ + 0{T^). 



(A43) 



At finite a and T ^ 1/a, the force becomes classical 

CB 



TF^E^{a), with 
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and 
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C(3) 



8a5/27r(a + 2Rf/^ 

In Eqs. (jA44p and (jA45p . we set R = 1, general expres- 
sions can be reconstructed by dimensional analysis. 
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